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1 Introduction 

The precise measurements in polarized M0ller scattering play a very impor- 
tant role in the modern polarimetry. The coincidence detection of the final 
electron pairs allows to essentially reduce a background of radiative effects [1] 
that accompany any charge particle scattering, including investigated reac- 
tion. However due to finite detector resolution it is impossible to remove all 
radiative event contributions out of the data. Moreover the additional virtual 
particle contributions can not be removed by any experimental cuts. Therefore 
to reach the appropriate accuracy we need to perform the radiative correction 
procedure consisting in the calculation of these effects within QED and esti- 
mate them numerically. 

The calculation of the lowest order QED radiative corrections (RC) to polar- 
ized M0ller scattering was already performed in [2,3] (see also the references 
therein). In [2] these corrections were calculated exactly (without the ultra- 
relativistic approximation) but unfortunately also without any experimental 
cuts. In [3] the similar corrections were performed for longitudinally polarized 
M0ller scattering within the ultrarelativistic approximation. The numerical 
analysis which is presented in [3] shows that RC to M0ller scattering are very 
sensitive to the missing mass (inelasticity) cuts. 

However, the consideration only the missing mass cuts during RC procedure 
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is insufficient. The realistic situation corresponds to taking into account the 
detector geometry that can essentially complicate the integration over the real 
photon phase space. In such situation an approach based on the Monte Carlo 
simulation of radiative events is the most adequate. The Monte Carlo gener- 
ators RADGEN [4] and ELRADGEN [5] for simulation of radiative events in 
deep inelastic and elastic lepton-nucleon scattering can serve as an examples 
of such approach. 

In present paper a new Monte Carlo generator MERADGEN 1.0 for the 
simulation of the radiative events in parity conserving part of polarized M0ller 
scattering is presented that can be used for the polarimetry purpose of present 
and future experiments, for example, in JLab and SLAG Naturally, as we are 
taking into consideration the parity conserving effects and the beam energy is 
rather small (from 1 to 45 GeV), the weak contributions (Z-boson exchange 
and so one) to this process are negligible. So, we restricted our calculation 
within QED theory only. 

The present paper is organized as following. In Section 2 the kinematics of the 
investigated process and the generation method are presented. The different 
contributions to the cross section that is responsible for real photon emission 
are considered in Section 3. Then in Section 4 the brief structure of the code 
is discussed, in Section 5 we described the input-output data. In Section 6 we 
explain how to run tests of our code, some conclusions are given in Section 7. 
Finally, in Appendices the 4-momenta reconstruction formulas, some lengthy 
formulas for RC and test output are presented. 



2 Kinematics and Method of Generation 



The lowest order contribution to the observable cross section of M0ller scat- 
tering reaction 



e(h, £ L ) + e(pi,r) L ) -> e(k 2 ) + e{p 2 ) (1) 

(in parentheses the 4-momenta and polarization vectors of electrons are pre- 
sented, and k\ = k\ = p\ = p\ = m 2 ) can be described by the standard set of 
Mandelstam variables: 



s = (ki +pi) 2 , t=(ki- k 2 ) 2 , u = (k 2 -pif 
s + t + u = Am 2 , 



(2) 



1 FORTRAN code MERADGEN 1.0 is available from |http:/ /www.hep.byTRC 
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while the beam and target (tjl) polarization vectors read: 



1 /s-2m 2 , o \ 
si = -7H «i - 2mpi), 



1 / s — 2m 2 \ 

where A s = s(s — Am 2 ). 

For the definition of the lowest order contribution to M0ller scattering it is 
enough to define an initial beam energy kio = (in Lab. system), a scatter- 
ing angle 9 C m (in CM system) of the detected electron with the 4-momentum 
k 2 and an azimuthal angle 0. The cross section does not depend on up to 
taking into account the detector geometry. Let us notice that E\ ah and cos 9qm 
can be expressed via the invariants in the following way: 



c — *2tti 2 

cos 6 CM = l + 2t/s = 1-22/, Eb ah = —z , (4) 

2m 

where y = —t/s. Therefore instead of the scattering angle and the beam energy 
for the definition of the Born cross section we can use s and t variables. 

For the radiative process with the real photon emission 



e(h, £ L ) + e(pi, r) L ) -> e(k 2 ) + e(p 2 ) + j(k) (5) 

(k 2 = 0) three new kinematic variables have to be defined. At first, following 
notations of ref. [3] we introduce an inelasticity v = A 2 — m 2 , where A = 
ki — k 2 + pi and A 2 is so-called missing mass square. Maximum value of the 
inelasticity 



st + Js(s - Am 2 )t(t - Am 2 ) 
v max = 1 s + t ( 6 ) 

can be defined from the kinematical restriction (see, for example, the Chew- 
Low diagram in [6]). This variable can be directly reconstructed from the 
data. To remove the contribution of the hard photon emission the events with 
v < v C ut are taking into account, where the value v cut is far less than v max . 
The second variable is defined by ti = (p 2 — Pi) 2 = (hi — k 2 — k) 2 . At last, 
the third variable should be choose as z = 2k 2 k. Notice that for the radiative 
process 

cosflcM = l + 2t/(s-v). (7) 
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To reconstruct the 4-momenta of all particles for radiative process in any 
system it is enough to determine of variables s, t, v, ti, z and the azimuthal 
angle <fi. As an example, in Appendix A the 4-momenta of electrons and real 
photon are expressed through these variables and presented in the center of 
the initial electron mass system. 

The simulation of the radiative events can be performed by the following 
algorithm: 

• For the fixed initial energy and t the non-radiative and radiative parts of 
the observable cross section are calculated. 

• The channel of scattering is simulated for the given event in accordance 
with partial contributions of these two (non-radiative and radiative) positive 
parts into the observable cross section. 

• The angle is simulated uniformly from to 2n. 

• For the radiative event the kinematic variables v, t% and z are simulated in 
accordance with their calculated distributions. 

• The 4-momenta of all final particles in the required system are calculated. 

• If the initial t has not a fixed value (i.e. simulated according to the Born 
probability distribution) then the cross sections have to be stored for 
reweighing. The t-distribution is simulated over the Born cross section, and 
realistic observed t-distribution is calculated as sum of weights, they are 
ratios of the observable and Born cross sections. 

Let us consider some important steps of simulation of the radiative events in 
more details. 



3 Non-radiative and radiative parts of the observable cross section 

Here we consider the observable cross section that has a form I" 2 !: 



where a is the Born contribution, a is an infrared divergency free sum of 
the contributions of the additional virtual particles and the "infrared" part 
of the real photon emission, crj? is the infrared divergency free part of the 
unobservable photon emission. The explicit expressions for each term of equa- 
tion (8) can be found in [3]. Since the expressions for the two virtual photon 
contribution as well as 5± and Sf in [3] contain the misprints (fortunately, 
it is not seriously reflecting on the numerical estimations presented in [3]), 



2 Here and later we consider the differential cross section a = da/dy 
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the additional virtual particle contributions together with 5± and <5f are also 
presented in Appendix B. 

Now we consider a R in more detail: before integration over the inelasticity and 
the real photon phase space it can be presented as 



^ v cut 

^(vcut) = ^ / dvj dF k (\M R \ 2 - 4F IR \M \ 2 ), (9) 
o 

where M R is a sum of matrix elements contributed to the real photon emission 
(see Appendix C for details), while M is a Born matrix element. The real 
photon phase space reads 

dTk = if^S((A - kf - m 2 ) = J- , dtldZ (10) 
where A(/c 1? k 2 ,pi, k) is the Gram determinant. At last, 



r R m 2 m 2 m 2 m 2 s — 2m 2 s — 2m 2 
-F IR = — + - T + — + - T + + + 

Z A Z{ V z V\ Z\V\ zv 

t — 2m 2 t — 2m 2 u — 2m 2 u — 2m 2 ,„ „ . 

+ + + — + — , 11 

Z\Z V\V ZiV V\Z 

where z\ = 2kk±, v\ = 2kp\ and u = s + t — 4m 2 . 

For separation of the cross section (8) into the radiative and non-radiative 
parts it is necessary to introduce a new fictitious parameter v m i n associated 
with missing mass square resolution. Then the equation (9) can be rewrite in 
the following way 



^ u min 

VF(vcut) = v r {v C ut,v mm ) + ^ J dvj dT k {\M R \ 2 - 4F IR \M \ 2 ) - 



2 ^cut 

J dvj dT k F IR \M \ 2 , (12) 



where 



O *-■(_; tit. ± *•> - 

ct f f f dz 

V r {Vcut,Vmin) = — / dv / dt ± / \M R \ 

167CS vL t L Z L ^-A(k l7 k 2 , Pl ,k) 
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dvd^dz y ' 



is the radiative part of the cross section. The limits of the integration over z 
are defined from the equation A(k±, k 2 ,pi, k) = 0, while t mm and t™ ax can be 
found from z mm = z max (see [6] for details). 

Since the observable cross section (8) is the sum of the radiative and non- 
radiative parts 



cr° hs {.v cut ) = a r (v cut , v min ) + o nr {v cuU v min ), (14) 
we immediately find that 

* nr (v cut ,v mm )=a° + a RV (v cut ) + ?L J dvj dT k (\M R \ 2 - AF IR \M \ 2 ) - 

o 

j dv f dT k F IR \M \ 2 . (15) 



Let us notice that the explicit formulae both for d 3 a r (v cut ,v m i n )/ dvdtidz in 
the equation (13) and <J nr {v cut ,v m i n ) allow us to start the generation of the 
radiative events. However to speed up the process of generation it is useful to 
perform the integration over z and t x analytically. So, the following analytical 
expressions 



do r {v cut, Vmin) _ f rf f ^ d 3 (J r (v cut , V min ) 

dv J J dvdtidz 

, d 3 a r ( 

Vcut i Vmin ) / -\ r\ 

ti J dvdtidz 



dvdU 



are incorporated in our Monte Carlo program. The analytical integration can 
be performed in a standard way (see, for example, [2,6] and references therein). 

At the end of this section it should be noted that according to the equation 
(14) the non-radiative and radiative contributions to the cross section depend 
on v cut and v m i n , while the observable cross section depends on v cut only. 
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program main 



urand I 




Fig. 1. The structure of the program MERADGEN 1.0 
4 The structure of the program and radiative event simulation 



The structure of the Monte Carlo generator MERADGEN 1.0 is presented in 
fig. 1. The main blocks mean: 

• merad_init - here we define all constants which are necessary for gener- 
ation; 

• grid_init — here we prepare the grids for generation of kinematic variables, 
really we approximate the theoretical curve by some sets of segments; 

• urand — random number generator (flat); 

• sig — the Born cross section and part of the virtual contribution (vertices 
and self energies); 

• xsbt — box contribution; 

• dcanc - the contribution with cancellation of the infrared divergency; 

• fsir - analytical cross sections da r (v cut ,v m i n )/dv, d 2 a r (v cu t,v m i n )/dvdti 
and d 3 a r (v cut , v min ) / dvdt x dz. 

For the convenience sake we split our programm into 4 FORTRAN-files: 

• run.f — main program for the event generation; 

• test.f — main program for the tests; 

• meradgen.f — the collection of main functions and subroutines of MER- 
ADGEN 1.0; 

• fsir.f — function fsir for calculation of the analytical cross sections that 
are presented above. 

For the events (tests) generation we need to run "make" ("make test") com- 
mand. 



Now let us consider the input-output data in more details. 
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5 Input-output data 



As an input data MERADGEN 1.0 uses 4- momentum of the virtual photon 
vpgen := k x — k 2 that generated in CM system externally, energy of electrons 
and degree of electrons polarization. There is only one variable itest in MER- 
ADGEN 1.0 that responsible for output. If itest := the output data are 
gathered in two common blocks (see file output. inc): 

common / variables / vgen , t lgen , zgen , weight , ich 

and 

common / vectors / vprad,phirad. 

Here vgen, tlgen and zgen are generated photonic variables v, t\ and z re- 
spectively, (they are necessary first of all for the test), weight is a ratio of 
the observable cross section to the Born one, variable ich shows the radiative 
(ich := 1) or non-radiative (ich := 0) scattering channel, the 4-momentum 
vprad := P2 — Pi and photonic 4-momentum phirad := k also defined in CM 
system. 

Here we have to do some remarks: 1) for non-radiative events vgen := 0, 
tlgen := t, z : = 0, vprad := vpgen and phirad := 0, and 2) as it was 
mentioned above, the variable v can be reconstructed experimentally and the 
events with hard photons v cu t < v < v max usually remove from the data. In 
order to speed up the process generation, we generated variable v from v m i n 
up to v cut = (s + t)/2 ~ v max /2. 

Now let us consider one sample of generation. It is well known that for the 
elastic process (1) the energy of the detected electron /c 2 o can be directly 
defined via the scattering angle 8cm and the initial electron beam energy in 
the following way: 



f£ab ^ CQS 2 l^itfLab (17) 

but when we deal with the real photon emission (5) the process becomes 
inelastic and, as a result, the "elastic" equation (17) is broken and we have 
the 2-dimensional distribution over and cos#cm- From fig. 2 one can see 
that this distribution has a sharp peak near an "elastic" line described by the 
equation (17). 
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Fig. 2. The 2-dimensional distributions of the scattering electron energy in lab. 
system and scattering angle 6 in CM system 

6 Test runs 



For our test runs we use the fact that if the events are simulated correctly 
their distributions over variables v, t\ and z must obey to the corresponding 
probability distributions: 



p{v) = 
P(ti) = 



1 da r (v cut ,v rn i n ) 
dv ' 
1 d 2 a r (v cut ,v m i n ) 



N v (J {Vcut-i^min) i 
d(7 (V cut , frran) 



Ni 



ti 



dvdt\ 



1 d 3 a r (v cut , 
N z dvdtidz 



ti 



dv 



d 2 a r (Vcut,Vmin) 

dvdti 



(18) 



Then for generation of p(v), p{t\) or p{z) distributions one has to put in the file 
test.f the value of variable itest such as: itest :— 1, itest := 2 or itest := 3, 
respectively, next to type "make test" and, at last, " ./test.exe" . The value 
rgen/rcalc, i.e. ratio of generated p-distribution to corresponding calculated 
cross section should be near unit. In Appendix D the test outputs for v,t±,z 
generation with P = 1 (see formula (B.2)), £ 6 Lab =45 GeV, 6 CM = 90°, 20 bins 
for the histogramming and 10 8 radiative events are presented. 

The simulated distributions of the photonic variables for the SLAC E158 ex- 
periment kinematic conditions and for the different degrees of polarization are 
presented in fig. 3 (all of parameters are noted there). We suppose v min rather 
small: v min = 2 x 10~ 2 E\ ah m. In the fig. 3 we can see clearly the divergent 
behavior of the distributions at v — > 0, z — > z mm ' max , t\ — > t™ m corresponding 
to the infrared singularity at the v m i n — > 0. Also it can be seen the physical, 
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Fig. 3. v, t±, z-histograms (points) and corresponding probability densities (curves) 
for E158 (SLAC) kinematic conditions (£^ ab = 45 GeV) 
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so-called z- and ^i-peaks of distribution p(t\). 

The other test with itest : = 4 consists in the cross-check of the accu- 
racy of vector reconstruction. So, using all 4-momenta we reproduce the 
value of generated invariant and compare them with the generated value. 
In this test we also calculate the photonic mass square as mlgamma : = 
phirad(4) 2 — phirad(l) 2 — phirad(2) 2 — phirad(3) 2 , which should be near zero. 
In last part of Appendix D we can see good coincidence reconstructed and 
generated invariants in different (random) kinematical points and very small 
values of photonic mass square. 



7 Conclusion 

In this paper the Monte Carlo generator MERADGEN 1.0 serving for the sim- 
ulation of radiative events in parity conserving longitudinally doubly-polarized 
M0ller scattering is presented. Following for the absolute necessity of both ac- 
curacy and quickness for our program we have developed the fast and high 
precise code using analytical integration wherever it was possible. MERAD- 
GEN 1.0 can be employed for the radiative corrections procedure in experi- 
ments with the complex detector geometry, such as SLAC E158 experiment 
and experiments of modern polarimetry (JLAB, SLAC). 
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Appendix A 4-momenta reconstruction 

The 4-momentum definition in the center of mass system of the initial electrons 
for M0ller process with real photon emission is shown in fig. 4 and can be 
presented in the form: 

h = (k w , 0, 0, |fci|), pi = (pio, 0, 0, -|pi|), 

k 2 = (k 20 ,k 21 ,k 22 ,k 23 ), P2 = (P20,P21,P22,P23), 

k = (k Q , h, k 2 , k 3 ), (A.l) 
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Fig. 4. Definition of momenta and angles in center-of-mass frame 



while their components can be expressed via the invariants and azimuthal 
angle (that is usually generated uniformly) in the following way: 



fcio = Pio = 7>Vs, N = \pi\ = 



ho 

P20 = 
P21 = 
P22 = 
P23 



S-V As • jl i v 7 ^ 

■7Tr> k21 = V^r cos< ?> k ™ = W — sm0, k 23 = 777^=, 



s — z 



2^' 



V A s AiA 8 sin + (4A 3 A 4 + sA 2 A 7 ) cos< 
4AiVA,A 3 
VAgAiAg cos0 — (4A 3 A 4 + sA 2 A 7 ) sin0 



4Ai V / AlA^ 



k 



fci = 



a/s(A 7 - A 2 A 4 ) 

2AiV / Aj 

f + 2; 



VAgAiAg sin0 + (4A 3 A 6 — sA 2 A 7 ) cosi 



4AiVATa; 



ko = 



VAsAiAs cos0 + (4A 3 A 6 - sA 2 A 7 ) sin0 
4AiV%A^ ' 
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y/s(X 7 + A 2 A 6 ) 
2AiV / Aj 



Here 



(A.2) 



Ai = (s - v) 2 - Asm 2 , \ 2 = 2t + s- v- 4m 2 , 

/ O \ O O / 0\/\ 

A3 = —st(s + t — v — Am ) — m v , A4 = s(s — v — Am ) — (s + 

A 5 = vz(s — v — z) — m 2 (y + z) 2 , A 6 = s(v — z) — v(v + z), 

A 7 = (s + 2ti - 2 - 4m 2 )Ai - A 2 A 4 , A 8 = I6A3A5 - X 2 . (A.3) 

As a result, the angles in fig. 4 can be expressed via the invariants in the 
following way: 



a s ^ 2 a s(A 7 + A 2 A 6 ) V^AiAs 

VA s Ai (w + ^jAiVAs sA 2 A 7 -4A 3 A 6 



Appendix B Additional virtual particle, 5± and 5f contributions 



The virtual contributions to M0ller scattering can be separated into three 
parts: 

a Y = a s + <r Vcr + <r Box , (B.l) 



where 1) <r s is a virtual photon self-energy contribution, 2) a Ver is a vertex 
function contribution, 3) <r Box is a box contribution. Now we consider each of 
them. 

(1) The contribution of the virtual photon self energies (including the photon 
vacuum polarization by hadrons) to the cross section looks like 



a Ana 2 

<T = 



Re(~±l(t) + U h (-t)) Ul + P)--(1-P)- 

v t ' L s u 



t 2 V t 

+(*<-> u). (B.2) 

Here P = PbPt, where Pb (Pt) is the beam (target) polarization, t) 
is the renormalized transverse part of the 7-self-energy [7] (this part in- 
cludes vacuum polarization by e, // and r charged leptons: in correspond- 
ing formula of [7] we should take a summing index / = e, /i,r). The 
hadronic part of the photonic vacuum polarization associated with light 
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quarks can be directly obtained from the data on process e + e — > hadrons 
via dispersion relations. Here we use parameterization of [8] 

ReU h (-t)^A + B\n(l + C\t\), (B.3) 

with updated parameters A,B,C in different energy regions. 
(2) For the contribution of the electron vertices we used the results of the 
paper [7] (see also references therein). We can obtain the vertex part as 



r Ver . 



2tt 3 



u 



s 2 i 



1 +P) (1 -P)— Ai(t,m 2 ) + (t <->u), 

s u J 



(B.4) 



where 



A ± (t, m 2 ) = -2 In §(ln 4 " l) + ^ 4 + ln ' 4 + ^ - 4 - ( R5 ) 

(3) Recalculated here expressions for the box cross section are slightly differ- 
ent from presented in [3] (we correct the misprints in the expression (16) 
of [3]) 



a Box = 



2a 2 



1 + P/2« 2 



■( In - — - In 

s v t \u\ 



u 



A 2 



^77)) 



\-P ( 2s\ s , \is 
^iX ln M ln A 2 



u 



(B.6) 



The expressions 5^ have the form: 



oL^ = l: — l s u - (r x + n z ) 



(77) 



°(77) 



2t 



l 2 sj + l X s-(ll + n 2 f-^- 



2 

p 

2 



(B.7) 



and logarithms look like 

S XL 

l s = In—, l x = In j. 



(B.8) 



It should be noted that vertex and box parts contain the infrared divergence 
through the appearance of the fictitious photon mass A. The infrared part 
from virtual cross section can be extracted in a simple way: 



*Y R = o v - a v (X 2 - a ) = In £(111 4- - (B.9) 
The correct expressions for 5f and 5f read 
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1_ 2 m { s(s + t-v ) 2 v (v -t) ) m 



1 j 2 ^max 



-\- 1 v max ^v ma x{v max £) 2 \t 

— ln 2 (l — Vrnax ) + In S - In ^ S — ^ S — - — Vmax ^ + 
t s -|- 1 v max t 2 



+ ln s — ^ „ ;= + ^ 1I ^ L + 2 

m s \t 



^max \ t ■ ( ^max 



Li, - Li 



±2 



t 



-Li 9 



s + t 



+L i 2 ™ - Li 



■"max 



2 



6 ' 



*f = - + (3 - 2/ r )/ m - (Z m - 1) ln - \ll - £ + 1. (B.10) 



Appendix C Matrix element of the real photon emission 



The sum of the matrix elements contributed to the real photon emission in 
M0ller process reads: 



m r = ]-u(k 2 )T m (k 2 , h)u(k 1 )u(p 2 )-f IJl u(p 1 ) + 
+ ju(k 2 )^u(k 1 )u(p 2 )r iJia (p2,pi)u(p 1 ) - 

- -u(p 2 )F m (p 2 , /ci)m(A:i)m(A;2)7^(Pi) - 

- — u(p 2 )^^u(ki)u(k 2 )T^ a (k 2 ,pi)u(pi). (C.l) 
The conjugate matrix element can be found in a simple way: 



M% = ^ L u(fci)r i/a (A;i, k 2 )u(k 2 )u(p 1 )~i v u(p 2 ) + 

1 _ _ 
+-«(A;i)7 1/ «(A;2)u(pi)r i/a (pi,p2)«(P2) - 

--^(^1)^(^1, p 2 )w(p2)«(pi)7^(^2) - 

— ^u(k{)^ v u{p2)u{p{)Y va {p u k 2 )u(k 2 ). (C.2) 

^2 

Here z 2 = z — s — t\ + 4m 2 and 



r, a (a,6) = (---) 7 ,- w -^ r , 
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7.. a 

By introducing 



r - (a ' h) = -(*& ~ y a n» --fir--™- (c - 3) 



S(a 1 , a 2 , a 3 , a 4 ) = Tr[aip(a 2 )a 3 A(a 4 )] 

S'(ai, a 2 , a 3 , a 4 , a 5 , a 6 , a 7 , a 8 ) = Tr[aip(a 2 )a 3 A(a 4 )a 5 p(a 6 )a 7 A(a 8 )], (C.4) 



where 



u(fci)w(A;i) = p(fci) = ^(fci + m)(l - Pb7 5 £l), 
u(pi)u(pi) = p(pi) = ^(p"i + m)(l - PtIsVl), 

u(k 2 )u(k 2 ) = A(fc 2 ) =k 2 + m, u(p 2 )u(p 2 ) = A(p 2 ), (C.5) 
the square of matrix elements reads: 



= -^S(T^ a (k 2 , fci), k ± , T ua (ki, k 2 ), A; 2 )S , (7 M ,j9i,7 1 ,,p 2 ) - 
h 

--S(r Ma (fc 2 , fci), fci, 7^ k 2 )S(^,p!, t va {pi,p 2 ),p 2 ) - 

-j-S(-f f ,,k l ,T ua (k 1 ,k 2 ),k 2 )S{T f , a (p 2 ,p 1 ),p 1 ,-f u ,p 2 ) - 
1 

; Pi 5 J- I'd (Pl,P2),P2) + 

+7^S'(iV(fc 2 , fci), fci,r„ a (fci,p 2 ),P2,7ji»Pi»7i'»fc2) + 
+7-^5'(r Ma (p 2 , fci), fci, f „ a (A;i, fc 2 ), fc 2 , 7m> Pi> 7^2) + 
+7— SfT^fca, fci), fci, 7„, p 2 , 7**»Pi» f „ a (pi, fo), fo) + 

+7 ^ in fa, 

(fci, fc 2 ), fc 2 , r MQ (fc 2 ,pi),pi,7 1 ,,p 2 ) + 
+— -S'(7 /1 ,fci,r i , a (fci,p 2 ) 

; P2; J- /iQ ip2,Pi),Pi,lu,k 2 ) + 

1 

+ — 5 ( (P2 , ^1 ) , ^1 , 7^ , ^2 , 7/* , Pi , T i/a (Pi , P2 ) , P2 ) + 
+ J~ S (7 M , fci , Iv, P2 , IV (p 2 , Pi ) , Pi , f ua (pi , fc 2 ) , fc 2 ) + 

+ 7 — 5 , (7^, fci, 7^ fa, r MQ (fc 2 ,pi) ; Pi 5 J- ua (P1,P2),P2) - 

r 1^2 
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\s(T^ a (p 2 , h), k 1 ,F va (k 1 ,p 2 ),P2)S('y li ,p 1 , 7„, k 2 ) 



v? 



1 Pi 1 J- VOL (pi,k 2 ),k 2 ) - 



uz 2 

Sij^ k 1 ,Tv a {ki,P2),P2)S{T iJia {k 2 ,p 1 ),p 1 ,~i v ,k 2 ) - 

uz 2 

— l - 2 S(-i^k l ^ u ,p 2 )SiT^ a (k 2 ,p 1 ),p u t ua (pi, k 2 ),k 2 ). (C.6) 

Z 2 



Appendix D Test output 



Here we present the results of the test as test.dat output file corresponding 
to: 

1) itest := 1 - the generation of p(v) distribution and comparison it with the 
analytical cross section corresponding to the first formula in (18) (here and 
below all of invariants v,t\,z are in GeV 2 ) 

itest=l 

v generation 

rgen is generated probability 
rcalc is calculated probability 
Ebeam=45 . GeV 

theta=90.0 degrees in CM system 

P=pb*pt=1.00 beam polarization times target polarization 
number of bins 20 

number of radiative events 100000000 
initial random number 12 



bin 




V 


rgen 


rcalc 


rgen/rcalc 


1 


0, 


.6998E-03 


376.4 


375.0 


1.004 


2 


0. 


, 1268E-02 


205.8 


205.4 


1.002 


3 


0. 


. 1826E-02 


143.4 


143.2 


1.001 


4 


0. 


.2382E-02 


111.3 


111.2 


1.002 


5 


0. 


.2936E-02 


92.07 


91.91 


1.002 


6 


0. 


.3490E-02 


79.40 


79.29 


1.001 


7 


0, 


.4043E-02 


70.70 


70.59 


1.002 


8 


0. 


.4596E-02 


64.45 


64.41 


1.001 


9 


0. 


.5148E-02 


60.02 


59.96 


1.001 


10 


0. 


.5701E-02 


56.81 


56.78 


1.001 


11 


0, 


. 6254E-02 


54.52 


54.57 


0.9991 


12 


0. 


.6806E-02 


53.26 


53.15 


1.002 


13 


0. 


.7358E-02 


52.40 


52.38 


1.000 


14 


0. 


.7910E-02 


52.27 


52.19 


1.001 


15 


0, 


.8463E-02 


52.69 


52.54 


1.003 
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16 


0. 


9015E-02 


53. 


.38 


53. 


,40 


0.9996 


17 


0. 


9567E-02 


54. 


.94 


54. 


.78 


1.003 


18 


0. 


. 1012E-01 


56. 


.53 


56. 


.69 


0.9973 


19 


0. 


1067E-01 


59. 


35 


59. 


.17 


1.003 


20 


0. 


1122E-01 


62. 


.15 


62. 


.29 


0.9978 



2) itest := 2 - the generation of p{t\) distribution and comparison it with the 
analytical cross section corresponding to the second formula in (18) 

itest=2 

tl generation 

rgen is generated probability 
rcalc is calculated probability 
Ebeam=45 . GeV 

theta=90.0 degrees in CM system 

P=pb*pt=1.00 beam polarization times target polarization 
number of bins 20 

number of radiative events 100000000 
initial random number 12 
v= 0.5749E-02 



bin 




tl 


rgen 


rcalc 


rgen/rcalc 


1 


-0. 


.2852E-01 


235.5 


201. 1 


1. 171 


2 


-0. 


2638E-01 


183.0 


313.5 


0.5836 


3 


-0. 


.2533E-01 


30.60 


28.38 


1.078 


4 


-0. 


2377E-01 


11.23 


10.86 


1.034 


5 


-0. 


2228E-01 


8.905 


8.711 


1.022 


6 


-0. 


.2081E-01 


9.842 


9.674 


1.017 


7 


-0. 


1933E-01 


14.89 


14.58 


1.021 


8 


-0. 


, 1763E-01 


76.89 


73.02 


1.053 


9 


-0. 


1689E-01 


73.09 


74.37 


0.9827 


10 


-0. 


. 1518E-01 


11.62 


11.45 


1.015 


11 


-0. 


.1369E-01 


6.749 


6.655 


1.014 


12 


-0. 


. 1212E-01 


7.542 


6.696 


1.126 


13 


-0. 


1095E-01 


5.605 


5.257 


1.066 


14 


-0. 


9398E-02 


2.392 


2.340 


1.022 


15 


-0, 


7939E-02 


1.711 


1.700 


1.007 


16 


-0, 


.6487E-02 


1.385 


1.381 


1.003 


17 


-0. 


.5041E-02 


1.206 


1.205 


1.001 


18 


-0. 


3594E-02 


1.142 


1. 132 


1.009 


19 


-0. 


.2129E-02 


1.250 


1.213 


1.030 


20 


-0. 


. 1505E-03 


11.30 


7.410 


1.525 



3) itest := 3 - the generation of p(z) distribution and comparison it with the 
analytical cross section corresponding to the third formula in (18) 
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itest=3 

z generation 

rgen is generated probability 
rcalc is calculated probability 
Ebeam=45 . GeV 

theta=90.0 degrees in CM system 

P=pb*pt=1.00 beam polarization times target polarization 
number of bins 20 

number of radiative events 100000000 
initial random number 12 

v= 0.5749E-02 
tl=-0.1437E-01 



bin 




z 


rgen 


rcalc 


rgen/rcalc 


1 


0. 


.9087E-02 


640.0 


755.7 


0.- 


8469 


2 


0. 


, 1043E-01 


55.53 


56.31 


0. 


9862 


3 


0, 


. 1159E-01 


2C 


).80 


20.85 


0. 


9977 


4 


0, 


. 1274E-01 


1C 


).47 


10.36 


1 


.011 


5 


0, 


.1388E-01 


5. 


989 


6. 


002 


0. 


9977 


6 


0, 


, 1501E-01 


3. 


898 


3. 


872 


1 


.007 


7 


0. 


, 1615E-01 


2. 


727 


2. 


717 


1 


.004 


8 


0, 


. 1729E-01 


2. 


074 


2. 


062 


1 


.006 


9 


0, 


, 1841E-01 


1. 


695 


1. 


692 


1 


.002 


10 


0, 


, 1955E-01 


1. 


508 


1. 


493 


1 


.010 


11 


0. 


.2068E-01 


1. 


430 


1. 


416 


1 


.010 


12 


0, 


.2181E-01 


1. 


447 


1. 


434 


1 


.009 


13 


0, 


.2295E-01 


1. 


559 


1. 


545 


1 


.009 


14 


0, 


.2408E-01 


1. 


771 


1. 


764 


1 


.004 


15 


0, 


.2521E-01 


2. 


151 


2. 


135 


1 


.008 


16 


0, 


.2634E-01 


2. 


765 


2. 


760 


1 


.002 


17 


0, 


.2748E-01 


3. 


917 


3. 


891 


1 


.007 


18 


0, 


.2862E-01 


6. 


302 


6. 


267 


1 


.005 


19 


0. 


.2979E-01 


13.29 


13.26 


1 


.002 


20 


0, 


.3112E-01 


108.4 


120.0 


0. 


9027 



4) itest := 4 - the cross-check of the accuracy of the vector reconstruction for 
5 random radiative events 

itest=4 

variable reconstruction 
Ebeam=45 . GeV 

theta=90.0 degrees in CM system 

P=pb*pt=1.00 beam polarization times target polarization 
number of radiative events 5 
initial random number 12 
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event= 1 

test v reconstruction 
v=0 . 612136E-02 reconstructed v from 4-vectors 
v=0.612136E-02 generated v 

test tl reconstruction 
tl=-0. 170339E-01 reconstructed tl from 4-vectors 
tl=-0. 170339E-01 generated tl 

test z reconstruction 
z=0 . 155138E-01 reconstructed z from 4-vectors 
z=0. 155138E-01 generated z 

m2gamma=-0 . 150213E-09 real photon mass square 



event= 2 

test v reconstruction 
v=0 . 861313E-03 reconstructed v from 4-vectors 
v=0.861312E-03 generated v 

test tl reconstruction 
tl=-0.221464E-01 reconstructed tl from 4-vectors 
tl=-0.221464E-01 generated tl 

test z reconstruction 
z=0 . 877074E-03 reconstructed z from 4-vectors 
z=0.877074E-03 generated z 

m2gamma=-0 . 231178E-11 real photon mass square 



event= 3 

test v reconstruction 
v=0 . 659541E-02 reconstructed v from 4-vectors 
v=0.659541E-02 generated v 

test tl reconstruction 
tl=-0. 162121E-01 reconstructed tl from 4-vectors 
tl=-0. 162121E-01 generated tl 

test z reconstruction 
z=0 . 689029E-02 reconstructed z from 4-vectors 
z=0.689029E-02 generated z 
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m2gamma=-0 . 71091 1E-10 real photon mass square 



event= 4 

test v reconstruction 
v=0 . 850422E-03 reconstructed v from 4-vectors 
v=0.850423E-03 generated v 

test tl reconstruction 
tl=-0 . 229561E-01 reconstructed tl from 4-vectors 
tl=-0.229561E-01 generated tl 

test z reconstruction 
z=0 . 196149E-03 reconstructed z from 4-vectors 
z=0. 196149E-03 generated z 

m2gamma= 0.105991E-12 real photon mass square 



event= 5 

test v reconstruction 
v=0 . 267845E-02 reconstructed v from 4-vectors 
v=0 . 267845E-02 generated v 

test tl reconstruction 
tl=-0 . 255805E-01 reconstructed tl from 4-vectors 
tl=-0.255805E-01 generated tl 

test z reconstruction 
z=0 . 154160E-02 reconstructed z from 4-vectors 
z=0. 154160E-02 generated z 

m2gamma=-0 . 316828E-11 real photon mass square 
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